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ABSTRACT 

We report measurements of the electron temperature at about a dozen 
locations in the extended emission-line regions of five active (Seyfert and radio) 
galaxies. Temperatures (^[oin] and TWjn) have been determined from both the 
I([OIII]A4363)/I([OIII]A5007) and I([NII]A5755)/I[NII]A6583) ratios. T [OIII] lies in 
the range (1.0 - 1.7) x 10 4 K. We find a strong trend for T[oin] to be higher than 
T[jv//], with the difference typically being « 5,000K. Because the critical density for 
colhsional de-excitation of the 1 T>2 level in Nil is lower than that of the same level in 
OIII, the deviations of the measured intensity ratios from those expected for T^qiii] 
= Tjjv//] in the low density limit are unlikely to result from colhsional de-excitation. 
The measured values of ^{oiii} and the differences between T[om] and Tyvr/i are 
very similar to those found in Galactic planetary nebulae. It is argued that the 
dominant form of energy input to the clouds is photoionization, but detailed modelling 
indicates that the temperature difference is too large to be accounted for in terms of 
photoionization of ionization-bounded clouds. We propose instead that both matter- 
and ionization-bounded clouds are present in the extended emission-line regions, with 
most of the [OIII] emission originating from a hot zone in the matter-bounded clouds 
and essentially all of the [Nil] from the ionization-bounded clouds. 
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1. Introduction 

High ionization, narrow emission lines are a denning characteristic of the nuclei of Seyfert 
galaxies and are often seen in radio galaxies. The gas emitting these lines is generally believed 
to be photoionized, though the origin of the ionizing photons is still debated. Spectroscopic 
studies have focussed mostly on the spatially unresolved nuclei, the narrow line emission of which 
originates from gas with a wide range of densities (10 2 < n e < 10 7 cm -3 ). This large spread of 
densities, together with the unknown distances of the individual gas clouds from the source of 
ionizing photons, leads to serious ambiguities in photoionization models. For example, collisional 
de-excitation of the X D2 level of OIII becomes significant for n e > 10 5 cm -3 , rendering the intensity 
ratio R[oni] = I(A4363)/I(A5007) an unreliable measure of electron temperature (Osterbrock 
1989). 

The high ionization gas is, however, often spatially extended in ground-based observations 
of both Seyfert (e.g. Mulchaey, Wilson & Tsvetanov 1996) and radio (e.g. Baum & Heckman 
1989) galaxies. Study of such extended gas provides two advantages. First, the density is low (n e 
< 10 3 cm -3 - e.g. Morganti et al. 1991; Tadhunter et al. 1994), so that collisional deexcitation 
is unimportant for most forbidden lines, the exceptions providing density diagnostics. Second, 
the geometric properties of the gas, including its distance from the nucleus, can be measured 
directly. These extended emission-line regions (EELRs) thus provide much simpler physical 
situations than the spatially unresolved narrow or broad line regions, and their excitation should 
be correspondingly easier to understand. 

Tadhunter, Robinson & Morganti (1989) have presented measurements of R[o/77] f° r 
twelve locations in EELRs around six active galaxies. These measurements imply electron 
temperatures in the range 12,800 < ^[oin] < 22,000K, whereas photoionization models which 
successfully account for most other line ratios predict T^qiii] < 11,000K. These high measured 
temperatures were confirmed by Storchi-Bergmann et al. (1996, hereafter SBWMB), who found 
10,000 < T^oin] < 17,000K for EELRs around five active (both Seyfert and radio) galaxies. 
Tadhunter, Robinson & Morganti (1989) concluded that either the EELRs have a heating source 
in addition to photoionization - plausibly cosmic rays or shocks - or the metal abundances 
are lower than the solar values assumed in the models. Binette, Wilson & Storchi-Bergmann 
(1996, hereafter BWSB) favored an alternative picture in which there are two populations 
of ionized clouds - a matter-bounded (MB) component responsible for most of the Hell and 
high-ionization forbidden line emission and an ionization-bounded (IB) component emitting low- 
to intermediate-ionization forbidden lines. In this model, the IB clouds are illuminated by an 
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ionizing continuum modified by absorption in the MB clouds. A new sequence of photoionization 
models was thereby obtained by varying the ratio, A M /i, of the solid angle subtended by the 
MB clouds at the ionizing source to that subtended by the IB clouds. The main success of 
the Am/ i sequence of models is that it provides a natural explanation for observed correlations 
between both the I([NeV]A3426)/I([OII]AA3727) and I([OIII]A5007)/I([OII]AA3727) ratios and the 
I(HeIIA4686)/I(H/3) ratio (BWSB). Such a correlation between the gaseous ionization, measured 
from the forbidden lines, and the I(HeIIA4686)/I(H/3) ratio cannot be understood in terms of the 
standard U (ionization parameter) sequence, since the HeIIA4686/H/3 ratio for IB clouds depends 
primarily on the slope of the ionizing continuum between 13.6 and 54.4 eV, and is relatively 
insensitive to U. Further, the electron temperature T[oni] is correctly predicted by the A M /i 
sequence, provided the thickness and ionization parameter of the MB component are appropriately 
selected. Finally, much stronger high ionization lines (e.g. [NeV]A3426 and CIVAA1549) are 
expected in the A M /i sequence than in the U sequence, again in accord with observations. 

In the A M /i sequence, the [OIII]A5007 and A4363 emission originates mostly in the MB 
clouds, while essentially all of [NII]A6583 and A5755 emission comes from the IB clouds. Because 
the MB clouds are hotter than the IB ones (see Figs 1 and 2 of BWSB), the temperature T\onn 
inferred from Rro/j/i should be considerably larger than the temperature T^//] inferred from 
R [NII] = I(A5755)/I(A6583). BWSB predict that T [om] should be ~ 5,000K larger than T [NII] in 
the A M /i sequence, but only slightly higher (by < 1,000K) in the [/-sequence. The temperature 
T[7V77] is predicted to be 9,200K for solar abundances in the A M /i sequence, given the value of U 
assumed by BWSB for the IB clouds. These predictions of the Am/i sequence models are in good 
agreement with temperature measurements of the NW cloud of Cygnus A, which give T^qiii] = 
15,000 ± 1,000K and T [NII] = 10,000 ± 600K (Tadhunter, Metz & Robinson 1994). 

In a typical Seyfert galaxy, I([OIII]A5007)/I(H/3) ~ 10, I(Ha)/I(H/3) ~ 3 and 
I([NII]A6583)/I(Ha) ~ 1. For T e = 10 4 K, R [NII] = 0.016 and R [on i] = 0.0064 (Osterbrock 1989). 
Thus, if T [NII] = T [OIII] = 10 4 K, [NII]A5755 and [OIII]A4363 should have roughly the same 
strength (I([OIII]A4363) ~ 1.3 x I([NH]A5755) and be about equally easy to detect with modern 
CCD spectrographs. However, this expectation is not borne out by our observations. 

The purpose of the present letter is to provide further measurements of T[nii] i n EELRs and 
evaluate the implications for the gaseous excitation. We find a strong tendency for Tr 0m i to 
exceed T[ NII ] by several thousand degrees and note that a similar trend, with a similar range of 
temperatures, is found in Galactic planetary nebulae. It is argued that the most likely explanation 
is that the EELRs contain both matter- and ionization-bounded clouds. 

2. Results 

In their study of nuclear and EELR spectra of 5 active galaxies, SBWMB report a detection 
of [NII]A5755 at only one location - the nucleus of Mkn 573. We have, therefore, measured upper 
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limits to the flux of this line at the 12 other locations where SBWMB obtained spectra. Two 
methods were used. In the first method, a spectrally unresolved artificial line was superimposed 
on the continuum and its flux adjusted by eye to provide an upper limit to the actual line flux. In 
the second method, the formal r.m.s. noise in the continuum was calculated and the 3 x r.m.s. 
upper limit on the flux of an unresolved spectral line was determined. In some cases, the upper 
limits for these two approaches agreed quite well, while in others the first method gave a higher 
value. We decided to adopt the more conservative upper limits from the first method. 

One concern was that the stellar template spectrum (which was subtracted from the observed 
spectrum prior to measuring the emission lines - see SBWMB) might, for some unknown reason, 
be a poor match near A5755 rest wavelength. If the template were too high, a real emission 
line could have been mistakenly subtracted away. Examination of of the two templates used by 
SBWMB showed them to be in excellent agreement in the spectral region in question and that 
A5755 corresponds to a local minimum. Since the templates match the stellar continuum very well 
almost everywhere else, we have no reason to doubt the validity of the template subtraction. 

The resulting upper limits to [NII]A5755, and the fluxes of [NII]A6583, [OIII]A4363 and 
[OIII]A5007, were corrected for obscuration, when present, using SBWMB's values of Ay, which 
were obtained assuming an intrinsic emission- line Balmer decrement of Ha/H/3 = 3.1. Correction 
for obscuration was not possible for the nucleus of ESO 362-G8, because a reliable measurement 
of H/3 could not be obtained due to the A type stellar continuum; the nucleus of this galaxy is 
omitted from the following discussion. The measurements of or limits to R[at//] and R[oni] were 
then converted to T[ NII ] and T^qiii] using equations 5.5 and 5.4 of Osterbrock (1989), respectively, 
assuming the low density limit. 

The results are shown in the left hand panel of Fig. 1, in which T^qiii] is plotted against 
T[nii]- It is immediately apparent that ^[oni] tends to be greater than T^//] by an amount 
which ranges up to rj 7,000K, but is more typically « 5,000K. This trend is found for the nucleus 
of NGC 526A and for most of the off-nuclear locations (only upper limits to both temperatures 
are available for the nuclei of PKS 0349-278 and PKS 0624-206). It is unlikely that this trend is 
related to collisional de-excitation of the l T>2 levels, since a) the critical densities of 8.6 x 10 4 
cm -3 (for Nil) and 7.0 x 10 5 cm -3 (for OIII) are much higher than expected for off-nuclear gas, 
and b) such collisional de-excitation effects would tend to depress [NII]A6583 at a lower density 
than [OIII]A5007, thus increasing R[nu] and the calculated value of T[ NII j relative to R[om] an d 
T [0///]- 

The right panel of Fig. 1 shows a similar plot for Galactic planetary nebulae using the 
measurements of Kingsburgh & Barlow (1994). The range of [OIII] temperatures is similar to that 
seen in the AGN and there is also a trend for T^^oni] t° exceed Tr^jn (see also Torres-Peimbert 
k, Peimbert 1977; Kaler 1986). The trends between T^yoni] an d ^[nii] i n planetary nebulae are 
generally considered to be excitation effects in a photoionized gas (Kaler 1986). The similar 
behavior of the AGN EELRs then strongly suggests that the same may be true for them. 
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3. Discussion 

The photoionization code MAPPINGS IC (Binette et al. 1993; Ferruit et al. 1997) has been 
used to calculate the expected values of Rfom] an d Rr/vJil i n various situations. We began by 
considering IB, isobaric, dust-free (/i = 0)0 clouds of solar metallicity (Z = 1) and low density (n 
= 1,000 cm -3 at the irradiated face) photoionized by continua of various shapes - blackbodies, 
power laws and a spectrum with a uv bump. A grid of models was generated by varying the 
ionization parameter U. The results are plotted in Fig. 2, in which the models represented by the 
long dashed (power-law ionizing spectrum) and continuous (blackbody-shaped ionizing spectrum) 
lines have [OIII]A5007/H/3 < 16, and are thus consistent with the observed values of this ratio, 
while the dotted extensions have [OIII]A5007/H/3 > 16, and are thus unacceptable. The models 
with [OIII]A5007/H/3 < 16 do show deviations from the Ttoni] = ^[nii] line, but these differences 
between T^yoiii] an d ^[nii] are much smaller than are observed. A power-law model with very low 
metallicity (Z = 0.2) provides [OIII] temperatures in general agreement with those observed, even 
without dust ({j, = 0), but again the difference between Tiqiii] and Tnvu] is too small. Also, such 
a low metallicity is implausible for the circumnuclear emission line regions discussed here. A high 
density (n = 60,000 cm -3 ) model lies to the right of the ^[oni] = ^[nii] nn e, confirming that 
collisional de-excitation tends to increase R[tv/7] relative to ~R[oni], which is opposite to the trend 
observed. 

As a second step, we have considered MB clouds. We adopt the model of BWSB, which 
contains two cloud populations - high ionization MB clouds and low ionization IB clouds. The 
two vertical dashed lines in Fig. 2 show the predicted behavior of R[om] an d R[7v/7] as a function 
of the parameter Am/i (defined in Section 1). The model marked 'A^/j a = -1.3' is the same as 
described in BWSB, and represents a varying proportion of MB clouds (with Umb = 0.04) and IB 
clouds (with Uib = 5.2 x 10~ 4 ). The model marked ''Am/i a = -1.1' has a harder continuum, 
Umb = 0.03 and \Jib = 1-8 x 10~ 4 . The predictions of these models are in agreement with the 
measured values in the nucleus of Mkn 573 and the NW cloud of Cygnus A and are consistent 
with most of the other points. These models were also shown by BWSB to be in agreement with 
the other optical line ratios. It should be emphasised that our models assume a single ionization 
parameter and total column for each of the matter- and ionization-bounded populations of clouds, 
which represents a gross oversimplification of a probably complex situation involving a continuous 
range of Umb, Ujb and the column density through the clouds. The model predictions should 
thus be considered as illustrative rather than definitive. Nevertheless, the key prediction (BWSB) 
of the matter- plus ionization-bounded model - that T[oiii] should exceed Tnvu] by ~ 5,000K - 
seems to be borne out by the data. 

Potential sources of uncertainty are non-collisional contributions to [OIII]A4363. Such 
contributors include charge transfer (e.g. Dalgarno &: Sternberg 1982) and recombination (e.g. 



4 /i is the dust-to-gas ratio of the cloud in units of the solar neighborhood dust-to-gas ratio. 
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Rubin 1986). The relative contribution of these effects is lower at higher temperatures because 
of the exponential temperature dependence of collisional excitation. Although no detailed 
calculations are available in the AGN context, and population of the 1 So level of OIII by these 
processes is not included in MAPPINGS IC, we infer from calculations with hot stellar ionizing 
sources (Kingdon & Ferland 1995) that the resulting error is less than 300K when T\qni] > 
10,000K. Therefore, recombination processes make a very minor contribution to [OIII]A4363 in 
the AGN context. 

Lastly, there is the question of whether sources of energy other than photoionization may 
contribute to the temperature difference between [OIII] and [Nil]. Such sources include shocks, 
which may be generated by outflowing winds or jets, and relativistic particles. While temperature 
stratification is expected in classical shock models (e.g. Binette, Dopita & Tuohy 1985), it is 
important to emphasise that the line ratios seen in the extended emission line regions under study 
are characteristic of photoionization by a hard spectrum (BWSB). Shocks with photoionized 
precursors can generate spectra which resemble those of Seyfert galaxies (Dopita & Sutherland 
1995), but then the temperature is expected to be similar to that found in classical photoionization 
models of IB clouds (see Fig. 7 of Dopita & Sutherland 1995). Relativistic particles may contribute 
heating (Ferland & Mushotzky 1984), but most of the off-nuclear regions we have studied do not 
have strong radio emission. Also, it is unclear why the relativistic particles would heat the OIII 
region preferentially over the NIL 

In summary, the range of electron temperature and the difference between T^qiii] and T^//] 
are very similar in planetary nebulae and the EELRs of active galaxies. These similarities do 
not favor a major contribution to the temperature from shock heating in the galaxies, despite 
the higher gas velocities there. More generally, the similarities between the two classes of object 
support the idea that the difference between the [OIII] and [Nil] temperatures is primarily related 
to the structure of the ionized clouds, rather than the presence of extra heating sources or a 
property of the ionizing source. Indeed, the presence of matter-bounded clouds in some planetary 
nebulae is demonstrated by differences between the two Zanstra temperatures derived from the 
nebular H I and He II lines (Osterbrock 1989, p 150), and also possibly by discrepancies between 
the Zanstra temperatures and the spectroscopic effective temperature of the central star (e.g. 
Kudritzki & Mendez 1989). 

ASW thanks J. P. Harrington for several valuable discussions. This research was supported 
in part by NASA under grants NAGW 4700 and NAG 81027 and by the Space Telescope Science 
Institute through grants G05411 and GO6006. LB is grateful to CNPq and CAPES of Brazil for 
financial assistance. 
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CAPTIONS TO FIGURES 

Figure 1 — Left: A plot of [OIII] temperature versus [Nil] temperature (in K) for nuclear and 
extended gas in active galaxies. Temperatures are plotted for all locations at which spectra were 
obtained by Storchi-Bergmann et al. (1996), except for the nucleus of ESO 362-G8 (see text) 
and 10" NW of the nucleus of PKS 0634-206, at which only a poor upper limit on T^//] could 
be obtained. The error bars are dominated by the uncertainties in the fluxes of [OIII]A4363 and 
[NII]A5755. The straight line represents Tfoni] = ^[nii]- The temperature sensitive [NIIJA5755 
line is detected in only the nucleus of Mkn 573, so the [Nil] temperatures are upper limits at all 
other positions. Right: A similar plot for planetary nebulae, constructed from the temperatures 
and their errors listed by Kingsburgh & Barlow (1994). It is notable that the range of T[oni] is 
similar in the two types of object and that there is a trend for T[ 0/// ] > T^ NII j in both types. 

Figure 2 — A plot of R [0 ///] = I(A4363)/I(A5007) versus R [NII] = I(A5755)/I(A6583). The 
diamonds represent the measured ratios in the nuclei of Mkn 573 and the NW cloud of Cygnus A 
(Tadhunter, Metz & Robinson 1994), while the error bars and upper limits represent our other 
measurements (see Fig. 1). The straight dot-dashed line is the locus ^[oni] = ^[nii] m the low 
density limit; the three asterisks represent temperatures of 8,500K, 10,000K and 15,000K (from 
lower left to upper right). All other lines represent model predictions calculated with MAPPINGS 
IC. Unless otherwise noted, the models are for ionization-bounded, dust-free, isobaric clouds of 
solar abundance and density 1,000 cm~ 3 . The two continuous lines are sequences of ionization 
parameter with a black-body shaped ionizing continuum having the temperatures indicated. For 
these models represented by the continuous lines, the ratio [OIII]A5007/H/3 < 16, in accord with 
observations. The dotted extensions to the continuous lines have [OIII]A5007/H/3 > 16. The 
four long dashed lines are sequences of ionization parameter for ionization-bounded clouds, with 
indicated abundance (Z) and dust content (/x), ionized by a power-law shaped ionizing continuum 
of index a = -1.3. These four models have [OIII]A5007/H/3 < 16 where the line is long dashes, 
but [OIII]A5007/H/3 > 16 for the dotted extensions. The long dashed curve to the right is for high 
density (60,000 cm~ 3 ) clouds. The long dash - short dash line towards the bottom is a sequence of 
ionization parameter for a broken power law ionizing spectrum (simulating a uv bump), following 
Mathews & Ferland (1987). Lastly, the two vertical short dashed lines to the left are models of 
the A M /j sequence with Z = 1, [i = 0.015 and power-law ionizing continua of index a = -1.3 and 
a = -1.1. The parameter A M n increases from 0.04 at the bottom to 16 at the top. 
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